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ABSTRACT 

We present deep Sparse Aperture Masking (SAM) observations obtained with the 
ESO Very Large Telescope of the pre-transitional disk object FL Cha (SpT=K8, 
d=160 pc), the disk of which is known to have a wide optically thin gap separating 
optically thick inner and outer disk components. We find non-zero closure phases, 
indicating a significant flux asymmetry in the K#-band emission (e.g., a departure 
from a single point source detection). We also present radiative transfer modeling 
of the SED of the FL Cha system and find that the gap extends from 0.06l^oi AU 
to 8.3±1.3 AU. We demonstrate that the non-zero closure phases can be explained 



1 Institute for Astronomy, University of Hawaii at Manoa, Honolulu, HI 96822, USA 

2 LESIA, CNRS/UMR-8109, Observatoire de Paris, UPMC, Universite Paris Diderot, 5 place Jules 
Janssen, 92195, Meudon, France 

3 Depart amento de Fisica y Astronomia, Universidad de Valparaiso, Valparaiso, Chile 

4 Departamento de Astronomia, Universidad de Chile, Camino del Observatorio 1515, Las Condes, San- 
tiago, Chile 

5 Departamento de Astronomia y Astrofisica, Pontificia Universidad Catolica de Chile, 7820436 Macul, 
Santiago, Chile 

6 UMI-FCA, CNRS / INSU France (UMI 3386) 

7 School of Physics, University of Sydney, NSW 2006, Australia 

8 Department of Physics and Astronomy, Macquarie University, NSW 2109, Australia 
9 Sagan Fellow, lcieza@ifa.hawaii.edu 

10 Millenium Nucleus 'Protoplanetary Disks in ALMA Early Science", Universidad de Chile, Casilla 36-D, 
Santiago, Chile 

11 Hubble Fellow 



-2- 



almost equally well by starlight scattered off the inner edge of the outer disk or 
by a (sub)stellar companion. Single-epoch, single-wavelength SAM observations 
of transitional disks with large cavities that could become resolved should thus be 
interpreted with caution, taking the disk and its properties into consideration. In 
the context of a binary model, the signal is most consistent with a high-contrast 
(AK# ^4.8 mag) source at a ^40 mas (6 AU) projected separation. However, 
the flux ratio and separation parameters remain highly degenerate and a much 
brighter source (AK# ~1 mag) at 15 mas (2.4 AU) can also reproduce the signal. 
Second-epoch, multi-wavelength observations are needed to establish the nature 
of the SAM detection in FL Cha. 

Subject headings: circumstellar matter — protoplanetary disks — stars: individ- 
ual (FL Cha) — planetary systems — techniques: interferometric 

1. Introduction 

Primordial circumstellar disks with inner cavities or gaps are known as "transitional" 
disks and can be identified by their peculiar Spectral Energy Distributions (SEDs), which 
have reduced levels of near- and/or mid-IR excesses with respect to the vast majority of Clas- 
sical T Tauri star (CTTS) disks. Nevertheless, not all transitional disk SEDs look the same 
and it is clear that not all inner holes and gaps are produced in the same way. Grain-growth, 
photoevaporation, and dynamical interactions with (sub)stellar companions can all result 
in transitional disk SEDs (see Williams & Cieza, 2011, for a recent review). Distinguishing 
between the processes that could be responsible for the unusual SEDs of transitional disks 
requires considerable information: SED shapes, accretion rates, disk masses, and multiplicity 
information (Najita et al. 2007). 

FL Cha belongs to a subclass of transitional objects known as "pre-transitional" disks, 
which is characterized by a pronounced "dip" in the mid-IR SED (Espaillat et al. 2007). 
The SEDs of pre-transitional disks can be reproduced with models presenting a wide (>5- 
10 AU) optically thin gap separating an optically thick inner disk from an optically thick 
outer disk. The inner and the outer disk components of pre-transitional objects have already 
been resolved by long-baseline interferometry observations in the near-IR (e.g., Olofsson et 
al. 2011; Tatulli et al. 2011) and the submillimeter (e.g., Andrews et al. 2011), respectively, 
confirming the gapped structure. These gaps are not necessarily empty. Deep polarized 
intensity images reveal a population of /xm-sized grains within some of the gaps (Dong et al. 
2012) and the accretion rates onto the stars in pre-transitional systems suggest that the inner 
disks are continuously replenished with material from the outer disks. Since pre-transitional 
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disks tend to lack stellar companions (Pott et al. 2010; Kraus et al. 2011), their gaps are 
best explained by the dynamical interaction of unseen substellar or planetary-mass objects 
embedded in the disk (Dobson- Robinson & Salyk, 2012). Recent high-contrast observations 
using the Sparse Aperture Masking (SAM) interferometric technique on Keck and the VLT 
have already identified companion candidates to two other pre-transitional disks, T Cha 
(Huelamo et al. 2011) and LkCa 15 (Kraus k Ireland, 2012). 

FL Cha is a K8 CTTS (Luhman 2007) in the Chamaeleon I molecular cloud, located 
at 160 pc (Whittet et al. 1997). As part of a program aiming the direct detection of 
young planets in transitional disks (see also Schreiber et al. 2013), we have obtained deep 
VLT-SAM observations of FL Cha. Our SAM data show non-zero closure phases indicating a 
second source of near-IR emission besides the central star. We also present radiative transfer 
modeling of the FL Cha SED, in order to better constrain the size of the gap, and discuss 
the possible nature of the newly identified source. 

2. SAM observations and data analysis 
2.1. Observations 

The VLT-SAM observations of FL Cha were performed on March 6, 2012, using the 
a 7-hole" mask (Tuthill et al. 2010) on the NAOS-CONICA (NaCo) Adaptive Optics system 
(Lenzen et al. 2003). The mask at the pupil-plane blocks most of the light from target 
and resamples the primary mirror into a set of smaller sub-apertures that form a sparse 
interferometric array with 21 baselines. SAM observations allow for exquisite calibration 
of the point spread function of the stellar primary and the suppression of speckle noise by 
the application of interferometric analysis techniques, such as the measurement of closure 
phases (the sum of the phases around any 3 triangle of baselines). The SAM technique is 
sensitive to companions in the ^0.5-5 A/D separation range (corresponding to ^30-300 mas 
for K^-band observations in the VLT) and can reach a contrast limit of AK^7 mag at A/D 
(Kraus & Ireland, 2012). 

The observing sequence consisted of multiple "visits" (4 observations of 25 frames of 
10 s in K#-band) of FL Cha, alternating with observations of the stars FI Cha, FK Cha, 
and 2MASS Jl 1082577-7648315 used as calibrators. During the same observing run, we also 
observed the close binary system RX Jl 106.3-7721 (Lafreniere et al. 2008) in the L'-band 
to validate our observing strategy and data reduction technique. The observations were 
reduced using the Paris SAMP pipeline as described by Lacour et al. (2011). Both FL 
Cha and RX J 1106.3-7721 show non-zero closure phases, indicting a departure from a single 
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point-source detection (see Figure 1). 

2.2. Binary model 

In order to constrain the properties of the SAM detections, we first performed a Monte 
Carlo Markov Chain analysis using a simple binary model with 3 basic parameters: two 
positional parameters, either separation and position angle (PA) or ARA and Adec, and the 
magnitude difference (Am). The x 2 maps of such models are shown in Figure 1 for both FL 
Cha and RX Jl 106.3-7721. Our procedure is identical to the one used by Schreiber et al. 
(2013). We find that the binary parameters for the RX Jl 106.3-7721 system are very well 
constrained by our observations and robust to the choice of priors (see Figure 2). The flux 
asymmetry signal in the closure phases of FL Cha is much weaker and results in much larger 
uncertainties in the model-derived parameters. While the PA is relatively well constrained, 
Am and separation remain degenerate and highly dependent on the choice of prior distribu- 
tions. This is illustrated in Figure 2, where we show the posteriori probability distributions 
for two different sets of priors: 1) uniform distributions for ARA, Adec, and Am and 2) 
uniform distributions for PA, Am, and the logarithm of the separation. The data favors a 
Am ^-4.8 mag source, but there is a long tail in the probability distribution extending to 
lower Am values. The posteriori distribution of the separation could be bimodal or unimodal 
depending on the adopted prior. The degeneracy between the Am and separation param- 
eters is a known problem for small separations (< A/D) in aperture masking observations 
(Pravdo et al. 2006, Martinache et al. 2009) and is clearly seen in the joint distributions 
shown in Figure 2. The complex x 2 surfaces in the joint distributions are consistent with 
two families of solutions: a relative bright source (AK ^1-3 mag) at ^15 mas (2.4 AU) and 
a much fainter one (AK ~4-5 mag) at ^30-40 mas (5.0-6.5 AU). The relative probability of 
these two solutions strongly depends on the choice of prior and can not be unambiguously 
estimated with the available data. Future H-band observations should provide the addi- 
tional resolution needed to break the Am-separation degeneracy and solve the two-solution 
ambiguity in the context of the binary model (see Schreiber et al. 2013). 

3. Disk model 

Espaillat et al. (2011) successfully reproduced the optical to 38 /xm SED of FL Cha 
adopting a simple model consisting of two vertical walls: one inner wall at the dust sub- 
limation distance of 0.04 AU that is responsible for the near-IR excess and an outer wall 
at 15 AU that reproduces the observed mid-IR excess. Here we use a more physical model 
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in order to provide further constraints on the size and location of the gap in the disk. We 
include photometry data at longer wavelengths (70 and 870 /xm), which are sensitive to the 
properties of the outer disk. 

3.1. Spectral energy distribution 

We constructed the FL Cha SED from the sources listed in Table 1. We have also 
obtained 870 fim photometry with the Atacama Pathfinder Experiment (APEX) using the 
LABOCA camera (Siringo et al. 2009). The observations were executed on July 24, 2012 us- 
ing the "wobbler on-ofP mode and reduced with the standard bolometer array data analysis 
package BoA0. FL Cha was detected with a signal to noise ratio of 12. The SED of FL Cha 
is also well sampled by a Spitzer-IKS spectrum covering the 5.2 to 38 fim region (Astro- 
nomical Observation Request #12696320). The 0.44 to 870 /im SED is shown in Figure 3. 
The optical and near-IR wavelengths were corrected by extinction, adopting Ay = 3.14 mag 
(calculated from the Rc-Ic color excess) and the extinction relations listed in Cieza et al. 
(2007). The median SED of CTTSs from Furlan et al. (2006) is shown for comparison. The 
sharp "dip" in the FL Cha SED around 15 fim indicates the presence of a wide gap in the 
disk. 



3.2. Radiative transfer modeling 

We model the observed SED using the Monte Carlo radiative transfer code MCFOST 
(Pinte et al. 2006). We parameterize the structure of the FL Cha disk with two independent 
components, an inner and an outer disk. Each component is described by the following 
parameters: the inner and outer radii (R in and R ow t, respectively) and the index 7 for the 
surface density profile (S(r) = Eio(r/10 AU) 7 ). The scale height as a function of radius is 
given by H{r) = Hi (r/10 AU)^. The grain size distribution has the form dn(a) oc a^da, 
between the minimum (a m ^ n ) and maximum (a max ) grain sizes. The width of the gap in 
the disk is simply given by: R ou ter,in - Rinner,out- The parameters a min , a max , p, ijj, 7, and 
Router, out) and the disk inclination are not easily constrained by the available data. We 
thus fix them to more or less "standard" values, which are listed in Table 2: For the dust 
composition, we follow Espaillat et al. (2011) and adopt a 40%/60% mixture of amorphous 
and crystalline silicates. 



1 http://www.apex-telescope.org/bolometer/laboca/boa/ 
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We adopt the fitting procedure described by Mathews et al. (2012), which uses the 
Levenberg-Marquardt x 2 minimization algorithm to calculate the numerical gradients of the 
X 2 function and determine the next point in the parameter space to be sampled until the 
algorithm converges to a x 2 minimum. We ran the search algorithm 10 times using different 
starting values to better sample the parameter space. Each of the runs results in a set of best- 
fit parameters, the distribution of which can be used to calculate the mean and associated 
uncertainty. 

The results of the 10 runs are listed in Table 2 and the structure of our best-fit model 
is shown in Figure 3. We find that the near-IR excess is best reproduced by a very narrow 
ring at the dust sublimation radius, extending from 0.04 to 0.06 AU, in agreement with the 
modeling results by Espaillat et al (2011). We also find that inner disks wider than AU 
significantly overproduce the observed near-IR fluxes and that the dip characteristic of pre- 
transitional disk SEDs disappears for inner disks wider than ^3 AU (for the given R ou ter,in 
value, resulting in gaps narrower than ~5 AU). In our model, the inner edge of the outer 
disk is located at 8.3±1.3 AU from the star. This distance is a factor of ^2 smaller than that 
estimated by Espaillat et al. The source of the discrepancy is unclear, but since our models 
share the same stellar parameters, it is likely to be related to the different parametrizations 
used for the disk structure (two vertical walls versus a full 3-D disk model) and the different 
grain size distributions adopted. Resolved submillimeter images with ALMA should be able 
to directly measure the cavity size and settle the discrepancy. Our SED model also constrains 
the properties of the outer disk. In particular, we find the FL Cha disk is relatively "typical" , 
with a mass of 7.5 M JUP and a scale height of ^1 AU at a radius of 10 AU. 



4. Possible interpretations of the SAM detection 

4.1. Background contamination 

The probability P(0, m) for an unrelated source to be located within a certain angular 
distance 6 from a particular target is given by P(0, m) — 1 — e -7rp ( m ) e2 , where p(m) is the 
cumulative surface density of background sources down to a limiting magnitude m (Brandner 
et al. 2000). FL Cha is 9.11 mag in K#-band. Therefore, the peak in the AK^ probability 
distribution correspond to an apparent K^— band magnitude of ^13.9. Since there are 
7482 stars brighter than = 13.9 mag within a 1 deg radius of FL Cha in the 2MASS 
catalog, the probability of a background source at a <0.050" separation is of the order of 
5 x 10 -6 . Background contamination can thus be discarded as a likely explanation for our 
SAM detection. 
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4.2. A stellar companion 

In Section 12.21 we found the SAM data is consistent with both a relative bright source 
(AK^l-3 mag) at a projected separation of ^2.4 AU or a much fainter one (AK ^4-5 mag) 
at ^5.0-6.5 AU. Evolutionary tracks of young low-mass objects are very uncertain, but can 
be used to address the nature of the putative companion and try to distinguish between a 
low-mass star, a brown dwarf, or a planet. According to the evolutionary tracks by Siess et al. 
(2000), FL Cha is a 0.6 M star and the former solution corresponds to a stellar companion 
with a mass in the <0. 1-0.3 M range. In the recent radial velocity (RV) monitoring study 
of Chamaeleon I objects performed by Nguyen et al. (2012), FL Cha showed a constant RV 
(16.9±1.1 km/s), close to the typical values in the region (~15.3±2 km/s), over the 1 month 
baseline of the study. While the precision and time baseline of the measurements are clearly 
not enough to rule out most stellar binaries, the RV measurements disfavor solutions with 
stellar companions at ^2.4 AU. As a reference, a 0.1 M companion to a 0.6 M star with 
circular edge-on orbit and a 2.4 AU semi-major axis has a period of 4.46 years and a velocity 
amplitude of ±2.5 km/sec. 



4.3. A brown dwarf or a protoplanet 

We now consider the nature of the source if AK# is close to the ~4.8 mag peak shown in 
the probability distributions from Figure 2. This peak corresponds to an absolute magnitude 
of ^7.9 at 160 pc. This renders our source a factor of 3 brighter than the protoplanet 
candidate identified by Kraus & Ireland (2012) within the gap of the LkCa 15 disk, and 
the "hot start" models by Chabrier et al. (2000) assign it a mass of ^15-20 M JUP , for an 
age of 1 Myr. Taken at face value, this would place the object at the bottom of the brown 
dwarf mass function. However, since the source seems to be located inside the gap of an 
accreting transitional disk, the inner disk of which is mostly depleted, it is reasonable to 
expect material from the outer disk to flow across the gap onto the inner disk and then the 
star. Under such circumstances, the low-mass object inside the gap should accrete most of 
the material being transported across the outer disk (Lubow & D'Angelo, 2006), resulting in 
significant accretion luminosity. SAM observations in H- and L-band would provide near-IR 
colors and help establishing whether the FL Cha detection is consistent with a brown dwarf 
or an actively accreting protoplanet surrounded by a disk, as seems to be the case in the 
LkCa 15 system. 
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4.4. Thermal emission or starlight scattered off the disk 

Our SAM detection is inconsistent with direct thermal emission from the inner disk. 
Our radiative transfer model suggests that the inner disk is <1 mas in diameter and the 
thermal contribution from the outer disk to the observed K-band flux is <0.01%. However, 
starlight scattered off the inner edge of the outer disk can in principle produce the observed 
signal. If the disk is highly inclined, the brightness of the projected rim would be asymmetric 
and could result in non-zero closure phases. To test this hypothesis, we generated a 2.2 fim 
ray traced image of the best-fit MCFOST model described in Section [3] and calculated the 
closure phases of the resulting image as seen at different position angles and inclinations 
(Figure [3j). We find that a disk inclined by ^60 deg from face-on with a position angl^l 
of ^150 deg results in closure phases that fit the SAM data almost as well as the best-fit 
binary model described in Section f2T2l does (x 2 ^140 vs. ^130). The difference in the \ 2 
results is mostly due to the fact that the disk model produces slightly smaller phase values 
compared to both the binary model and the SAM data. However, we emphasize that the 
2.2 /im model image depends on disk properties that are poorly constrained, such as the 
grain size distribution, the dust composition, and the detailed structure of both the inner 
and the outer disk. These disks parameters could be varied to improve the fit further. 
Pending more constraints (i.e., resolved images), the current model suffices to show that 
a disk (which also fits the SED) can produce the necessary closure phase signal. This 
demonstrates that the circumstellar disk must be taken into consideration when interpreting 
SAM data of transitional disks with inner cavities that are large enough to become resolved. 
High-resolution submillimeter imaging with ALMA would reveal the exact orientation and 
size of the cavity in the FL Cha disk. Such imaging should establish whether our SAM 
detection is in fact consistent with starlight scattered off the disk or, on the contrary, if the 
source is located well inside the hole, as in the case of the LkCa 15 system. Measuring the 
K#-H and K#-L colors of the source would also help testing the scattered light hypothesis. 
Looking for orbital motion from multi-epoch SAM observations should provide the ultimate 
test to distinguish between a companion and the scattered light scenario. 



5. Summary and Conclusions 

From VLT-SAM observations, we have identified a near-IR flux asymmetry in the pre- 
transitional object FL Cha. By modeling its SED, we find that the gap in its disks extends 



2 The position angle is the direction defined by the intersection between the plane of the disk and the 
plane of the sky. 



-9- 



from 0.06to!oi AU to 8.3±1.3 AU. We have considered several potential possibilities for the 
nature of the source: a low-mass star, a brown dwarf, a protoplanet, thermal emission or 
starlight scattered off the disk, or a background object. Only direct thermal emission from 
the inner disk and background contamination can be ruled out. We find that light scattered 
off the inner edge of the outer disk can result in closure phases that fit the SAM data almost 
as well as binary models do. Single-epoch, single-filter SAM observations of transitional 
disks should thus be interpreted very cautiously, taking the disk into consideration. In the 
context of a binary model, the closure phases are most consistent with a AK# ^-4.8 mag 
source at a 6 AU projected separation, but a much brighter one (AK^ ~1 mag) at 2.4 
AU can also reproduce the signal. H- and L-band SAM observations of FL Cha are highly 
desirable to 1) test for orbital motion, 2) break the current degeneracy between Amag and 
separation in the binary model, and 3) provide near-IR colors to help distinguishing between 
the possible explanations. Resolved submillimeter images with ALMA are needed to establish 
the orientation of the system and provide a direct measurement of the inner cavity size to 
test the scattered light hypothesis and to better constrain the properties of the outer disk. 

L.A.C was supported by NASA through the Sagan Fellowship Program. M.R.S., A.J., 
S.C., and F.M. acknowledge support from the Millenium Science Initiative, Chilean Ministry 
of Economy, Nucleus P10-022-F. G.S.M. acknowledges NASA/JPL funding support through 
grant RSA-1369686. A.L.K. was supported by NASA through the Hubble Fellowship pro- 
gram. 
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Table 1. FL Cha photometry data 



Wavelength 


Flux 


Flux 


Error" 


Telescope 


Reference 6 


(//m) 


(mJy) 


(mag) 


(mJy) 






0.44 


1.74E-01 


18.44 


30% 


ground based 


1 


0.55 


1.16E+00 


16.28 


30% 


ground based 


1 


0.65 


3.19E+00 


14.96 


20% 


ground based 


1 


0.80 


1.10E+01 


13.41 


20% 


ground based 


1 


1.25 


5.45E+01 


11.73 


15% 


2MASS 


2 


1.66 


1.17E+02 


9.90 


15% 


2MASS 


2 


2.20 


1.51E+02 


9.11 


15% 


2MASS 


2 


3.6 


9.66E+01 


8.66 


10% 


Spitzer 


3 


4.5 


8.20E+01 


8.35 


10% 


Spitzer 


3 


5.8 


5.60E+01 


8.28 


10% 


Spitzer 


3 


8.0 


3.95E+01 


8.02 


10% 


Spitzer 


3 


24 


1.01E+02 


4.67 


10% 


Spitzer 


3 


70 


3.01E+02 




15% 


Spitzer 


3 


870 


3.00E+01 




15% 


APEX 


4 



Note. — a the optical and near-IR uncertainties are dominated by the extinction correc- 
tions. ^References: (1) Gauvin & Strom (1992); (2) Skrutskie et al. (2006); (3) Spitzer 's 
Gould Belt Catalog. The 3.6-8.0 fim photometry have already been published by Cieza et 
al. (2009). The 24 and 70 /xm data have not been published before; (4) this work. 
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Table 2. Disk and Stellar Parameters 



Parameter 


value 


error 


Stellar parameters 


Stellar T eff [K] 


3850 


fixed 


Stellar Luminosity [L ] 


0.4 


fixed 


Stellar mass [M ] 


0.6 


fixed 


Distance [pc] 


160 


fixed 


Inner and outer disk parameters 


Inclination [deg] 


60 


fixed 


Grain size distribution slope, p 


-3.5 


fixed 


a ram [/^ m ] 


0.005 


fixed 


a max [M m ] 


3900 


fixed 


Surface density exponent, 7 


-1 


fixed 


Flaring exponent, ip 


1.1 


fixed 


Inner disk parameters 


Scale height at 10 AU Hio,inner [AU] 


0.16 


0.05 




5xl0" 4 


5xl0" 4 


Rmner,m [AU] 


0.04 


0.01 


Rmner,owt [AU] 


0.06 


0.05 


Outer disk parameters 


Scale height at 10 AU Hi , o ^ er [AU] 


1.2 


0.7 


Mass diskiOUter [Mjup] 1 


7.5 


1.0 


^outer,in [AU] 


8.3 


1.3 


^outer,out [AU] 


100 


fixed 



Assumes a gas to dust mass ratio of 100. 
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Fig. 1. — Top Panels: Closure phase as a function of baseline from our VLT-NACO/SAM data 
indicating significant flux asymmetries in both FL Cha (left) and RX J1106. 3-7721 (right). Lower 
Panels: x 2 maps resulting from fitting a binary model to the closure phases. The 5-a contours 
indicate the location of the best-fit companions. 
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Fig. 2. — Top Rows: the posteriori distributions of the 3 parameters in the binary models 
(Amag, separation, and PA) for FL Cha (left, using two sets of priors as described in Section l2?2|) 
and RX J1106. 3-7721 (right). Bottom Row: The joint distributions of separation and Amag. For 
FL Cha, these two quantities remain degenerate and highly dependent on the choice of priors. 
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Fig. 3. — Top Panel: The extinction-corrected SED of FL Cha showing a pronounced "dip" at 
mid-IR wavelengths, the defining feature of pre-transitional disks (left) and the disk structure of 
our best-fit SED model (right). Middle Panel: Quality of fit to the SAM data for our best-fit disk 
model as a function of disk PA, compared to the best-fit binary model and a single point source 
(left). K^-band ray traced image of the best-fit disk model (right). Bottom Panel: UV coverage 
of the SAM observations on FL Cha indicating the sign and value of the phases (left). The phases 
resulting from the best-fit disk model adopting the same UV coverage as in the real data (right). 



